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LONG-TERM  GOALS 

Our  goal  is  to  provide  a  complete  simulation  code  that  will  represent  and  predict  the  sediment  transport 
and  bed  features  on  the  continental  shelf  at  user-specified  resolution  by  using  state-of-the-art 
algorithms  for  the  physics  and  numerics  of  the  simulation  code. 

OBJECTIVES 

Our  primary  objective  is  to  simulate  the  ripple  climate  on  the  bed  of  the  inner  shelf  at  depths  on  the 
order  of  20  m  and  over  domains  ranging  from  centimeters  to  kilometers  in  support  of  the  Ripples  DRI 
experiments  and  analyses.  Our  secondary  objective  is  to  use  simulation  to  better  understand  the 
physics  of  ripple  formation  and  sediment  transport  in  this  environment. 

APPROACH 

In  order  to  simulate  the  formation  and  evolution  of  benthic  ripples,  we  employ  a  large-eddy  simulation 
code  with  suspended  sediment  and  moving-bed  modules.  The  suspended-sediment  module  follows 
that  of  Zedler  and  Street  (2001,  2006),  and  is  incorporated  into  the  parallel  large-eddy  simulation  code 
PCUI  (Cui  and  Street,  2001;  2004)  in  order  to  perform  high-resolution  simulations  of  sediment-ripple 
dynamics  on  parallel  computers.  We  have  employed  several  modifications  to  the  implementation  of 
Zedler  and  Street,  most  notably  the  implementation  of  the  bottom  boundary  condition  for  sediment, 
which  does  not  require  knowledge  of  the  near-bed  turbulence  parameterization.  This  is  outlined  in 
Chou  and  Eringer  (2008a).  Our  moving  bed  module  responds  to  the  high-resolution  shear  stress 
imposed  by  the  fluid  at  the  bed  and  moves  the  bed  position  Zb  according  to  the  formula  (Gessler  et  ah, 
1999;  Wuet  ah,  2000) 
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where  p'  is  the  porosity  of  the  bed  and  q  is  the  bed-load  transport,  whieh  we  determine  through  an 
empirieal  formula  similar  to  that  of  Meyer-Peter  and  Muller  (1948),  and  depends  on  the  near-bed  shear 
stress.  Beeause  we  employ  terrain-following  eoordinates,  movement  of  the  bed  requires 
implementation  of  numerieal  algorithms  on  moving  grids  in  a  manner  similar  to  the  method  of  Hodges 
and  Street  (1999),  in  whieh  the  boundary-following  grid  is  updated  in  response  to  the  moving  free 
surfaee. 

WORK  COMPLETED 

We  have  developed  a  moving  eurvilinear  grid  formulation  that  enables  simulation  of  the  Navier-Stokes 
equations,  along  with  the  transport  equations  for  the  suspended  sediment,  on  a  grid  that  moves  in 
response  to  moving  boundaries.  The  method  guarantees  eonsistency  with  eontinuity,  in  that  if  the 
suspended  sediment  eoneentration  is  uniform,  it  remains  uniform  in  response  to  flow  and  transport  in 
the  absenee  of  erosion  or  settling. 

RESULTS 

We  outline  the  methodology  behind  the  moving  grid  formulation  and  present  results  that  demonstrate 
the  potential  ramifieations  of  employing  methods  that  are  not  eonsistent  with  eontinuity  (Gross  et  ah, 
2002).  Details  of  the  formulation  ean  be  found  in  Chou  and  Fringer  (2008b).  On  a  two-dimensional 
moving,  generalized  eurvilinear-eoordinate  finite-volume  grid,  as  shown  in  Figure  1,  the  equation 
governing  transport  of  a  scalar  field  C  in  grid  cell  ij  is  given  by 

j,{j„c,Mv-u.)c,-(u-uXc..+{v-v.)c. =0,  (2) 

where  e,  w,  n,  and  5  subscripts  denote  the  east,  west,  north,  and  south  faces  of  the  two-dimensional 
computational  cell  with  volume  Jy,  U  and  V  are  the  contravariant  volume  fluxes  (in  units  of  m^  s'^)  in 
the  transformed  curvilinear  coordinates,  and  the  g  subscript  denotes  the  contravariant  volume  flux 
associated  with  grid  movement.  The  time -rate  of  change  of  the  cell  volume  is  governed  by  the 
continuity  equation  in  the  moving  curvilinear  coordinate  system, 

incompressibility,  i.e. 

U  -U  +V  -V  =0,  t4) 
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equation  (3)  leads  to  the  equation  governing  discrete  conservation  of  space. 
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dJjj 

dt 


-u  +v 

g,n 


V  =0. 

g,s 


(5) 


2 


In  order  for  equation  (2)  to  be  eonsistent  with  eontinuity,  or  CWC,  substitution  of  a  uniform  sealar 
eoneentration  field  into  its  time-diserete  form  must  yield  the  time-diserete  form  of  (5).  In  general,  this 
is  not  a  straightforward  task  and  speeial  eare  must  be  taken  when  diseretizing  equations  (2)  and  (5)  in 
order  to  ensure  CWC.  A  simple  example  of  violation  of  CWC  ean  be  demonstrated  by  diseretizing 
equation  (2)  with  the  seeond-order  Adams-Bashforth  method,  sueh  that 
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and  diseretizing  equation  (5)  with  the  first-order  forward  Euler  method,  sueh  that 
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The  diseretization  (6)  is  not  eonsistent  with  (7)  beeause,  upon  substitution  of  the  eonstant  sealar  field 
C=Co  into  equation  (6)  and  assuming  eontinuity  of  fluid  (4),  one  obtains 
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whieh  violates  CWC  beeause  it  is  not  identieal  to  equation  (7). 

While  CWC  would  be  satisfied  if  both  equations  (2)  and  (5)  were  diseretized  with  the  forward  Euler 
method,  diseretization  (7)  violates  eonservation  of  spaee  unless  speeial  eare  is  taken  in  eomputing  the 
eontravariant  volume  fluxes  assoeiated  with  grid  movement.  This  is  demonstrated  in  Eigure  2,  whieh 
shows  that  in  one  dimension,  equation  (7)  would  yield  the  eorreet  new  volume  at  the  new  time  step. 
However,  the  two-dimensional  form  yields  a  value  for  the  new  volume  that  is  smaller  than  the  true 
value  as  would  be  eomputed  from  the  eell  vertiees  whieh  are  updated  independently  of  equation  (7). 
Chou  and  Eringer  (2008b)  developed  a  teehnique  to  avoid  this  problem  by  requiring  that  fluxes 
assoeiated  with  grid  movement  are  always  ealeulated  at  time  step  n+1/2,  rather  than  time  step  n. 

As  an  example  of  the  effects  of  not  ensuring  CWC  in  a  simulation  of  transport  over  a  moving  bottom, 
we  employed  a  simulation  with  a  bed  that  moves  according  to  a  pre-specified  function  (rather  than  due 
to  a  dynamic  bed  model).  With  an  initially  uniform  concentration  field  and  in  the  absence  of  settling 
or  erosion,  Eigure  3(a)  shows  that  a  CWC  method  guarantees  that  the  concentration  remains  uniform 
as  the  bed  oscillates  in  time.  However,  Eigure  3(b)  shows  that  a  non-CWC  method  leads  to 
exponential  growth  of  the  near-bed  concentration  field  because  the  non-CWC  method  effectively 
introduces  a  source  of  mass  in  the  near-bed  cells  whose  volume  changes  considerably  due  to  the  grid 
motion.  While  this  case  shows  pronounced  effects  of  violating  CWC,  the  effects  in  a  real  suspended 
sediment  simulation  are  not  as  evident  because  of  the  transient  nature  of  the  flow  in  response  to 
settling  and  erosion.  However,  as  shown  in  Eigure  4,  the  suspended  sediment  concentration  in  a  CWC 
simulation  differs  significantly  from  that  in  a  non-CWC  simulation,  particularly  in  the  near-bed 
suspended  sediment  concentration  field,  which  is  much  higher  for  the  non-CWC  case. 


3 


Figure  1:  Depiction  of  a  two-dimensional  finite  volume  grid  cell  in  a  generalized  curvilinear 

coordinate  system. 


Figure  2:  Example  of  grid  movement  when  space  is  conserved  in  one  dimension  (a)  and  when  space 
is  not  conserved  in  two  dimensions  (b)  using  the  discretization  (7).  Figure  (b)  shows  that  the  shaded 
area  incorrectly  underestimates  the  actual  area  of  the  new  cell. 
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Figure  3:  The  influence  of  a  moving  grid  on  a  uniform  concentration  field  for  a  CWC  method  (left 
panels)  and  a  non-CWC  method  (right  panels).  T  is  the  period  of  oscillation  of  the  bed,  and  the 

velocity  field  results  only  due  to  bed  movement 
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Figure  4:  Suspended  sediment  concentration  over  a  growing  ripple  field  in  a  wave-induced  current 
that  oscillates  with  period  T.  The  left  column  uses  a  CWC  method  while  the  right  column  does  not. 
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